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The development of the magneto-optical 
trap revolutionized the fields of atomic 
and quantum physics by providing a 
simple method for the rapid production 
of ultracold, trapped atoms. A sim- 
ilar technique for producing a diverse 
set of dense, ultracold diatomic molec- 
ular species will likewise transform the 
study of strongly interacting quantum sys- 
tems, precision measurement, and physi- 
cal chemistry. We demonstrate one- and 
two-dimensional transverse laser cooling 
and magneto-optical trapping of the polar 
molecule yttrium (II) oxide (YO). Using a 
quasicycling optical transition we observe 
transverse Doppler cooling of a YO molec- 
ular beam to a temperature of 5 mK, lim- 
ited by interaction time. With the addi- 
tion of an oscillating magnetic quadrupole 
field we demonstrate a transverse magneto- 
optical trap and achieve temperatures of 
2mK. 

Over the past quarter century, the magneto-optical 
trap (MOT) has been extended to two dozen atomic 
species (7). This abundance of species makes ultra- 
cold atomic systems a powerful tool for studying a wide 
range of phenomena, from quantum-degenerate gases, 
physics beyond the Standard Model (2), and strongly 



correlated systems (3), to applications in quantum in- 
formation (4) and simulation, quantum sensing, and ul- 
traprecise optical clocks (5). Ultracold polar molecules, 
with their additional internal degrees of freedom and 
complex interactions, yield even richer phenomena and 
open the door to new fields such as ultracold chem- 
istry (6). 

Recently, many techniques have been developed 
for producing cold and ultracold samples of polar 
molecules. Magneto-association and adiabatic trans- 
fer (7) of ultracold atoms can produce ultracold samples 
of polar molecules; though this technique is currently 
limited to bialkali species. Buffer gas cooling (8) and 
molecular beam slowing techniques, such as Stark (9) 
and Zeeman (70) deceleration, can produce molecules 
cold enough to load into conservative traps. Further 
cooling of these trapped samples to temperatures below 
10 mK via evaporative or sympathetic cooling remains 
technically challenging. Optical cooling has been pro- 
posed (77, 12) and recently realized (73). A MOT (72) 
would be the ideal tool for producing ultracold trapped 
samples of diatomic molecules, much as it is for atoms. 

A MOT gains its utility by combining a spatially de- 
pendent trapping force with a fast dissipative cooling 
force. With 1/e cooling rates on the order of 10 5 s _1 , 
warm atoms can be cooled, trapped, and compressed 
in a few milliseconds over length scales of less than 
1 cm. To achieve these fast cooling rates, more than 
10 4 optical photons must be scattered at rates of more 
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than 10 6 s -1 , requiring a highly closed electronic tran- 
sition. The additional vibrational and rotational degrees 
of freedom present in molecules makes creating a closed 
cycling transition difficult in these systems. Neverthe- 
less, for certain molecules it is possible to construct 
quasicycling transitions with only a minimal increase in 
laser complexity (11, 12). One-dimensional transverse 
laser cooling (73) and longitudinal slowing (14) of a 
SrF beam have been demonstrated. Opto-electric cool- 
ing of CH3F molecules has resulted in temperatures as 
low as 29 mK (75). Here we demonstrate a new type of 
MOT, suitable for quasicycling transitions available in 
molecules. We implement this system for the molecule 
yttrium (II) oxide (YO) and observe both viscous cool- 
ing and a magneto-optical spring force in a 2-D MOT. 

We begin by presenting the level structure of YO 
and describe the procedure for creating a quasicycling 
transition. YO has a single naturally abundant iso- 
topomer, 89 Y le O, and relatively simple hyperfine struc- 
ture with nuclear spins Iy = 1/2 and Iq = 0. The 
main cooling transition proceeds on A 2 £ — > A 2 Hi/ 2 
at 614 nm, as shown in Figure 1A. The A 2 Yli/ 2 has a 
radiative lifetime of -f^ 1 — 33 ns (76), allowing for 
fast optical cycling. Diagonal Franck-Condon factors 
limit the vibrational branching of A 2 Hi/ 2 (17). Only 
two additional lasers at 648 and 649 nm to repump 
the v" = 1 and v" = 2 vibrational levels are needed 
to limit vibrational branching loss to <10~ 6 . The 
\A, v'; J' = 1/2, F', +) manifold forms a highly closed 
transition with the \X, v"; N" = 1, G", F", -) mani- 
fold as shown in Fig. IB, due to parity and angular mo- 
mentum selection rules (78). The magnetic field depen- 
dence of the ground manifold is shown in Fig. 1C. Here, 
the ground electronic state is labeled by Hund's case bps 
(19) quantum numbers \X, v"; N" , G", F",p), with vi- 
brational level v", rotational level TV", intermediate 
quantum number G" formed by coupling of electron 
and nuclear spin G = S + I, and total angular moment 
F" (20). The excited A 2 U 1/2 state is labeled by Hund's 
case a quantum numbers, \A, v'; J', F',p), where J' is 
the total electronic angular momentum, and F = J + I. 
The lowest rotational state J' = 1/2 is split into two 
states of opposite parity (p = ±) separated by a A- 
doubling of 4.5 GHz (77). This cooling scheme pro- 
vides a rotationally closed transition at the expense of 
repumping of the G" = 0, 1 hyperfine ground states. 
The repumping of hyperfine levels within each vibra- 
tional level can be achieved with a single laser by creat- 
ing frequency shifted sidebands using an acousto-optic 



modulator (18). In order to maintain optimal photon 
scattering rates for laser cooling, we destabilize optical 
dark states that form in the ground state Zeeman man- 
ifolds (27) by modulating the polarization of the cool- 
ing light between er + and o~ with a voltage-controlled 
waveplate (Pockels cell). The modulation rate should be 
similar to the optical pumping rate, which in our case is 
on the order of several 10 6 s _1 . 

This laser configuration for creating a quasicycling 
transition also lends itself naturally to creating a MOT. 
Figure ID shows the simplified YO level structure in- 
volved in the MOT. In contrast to a typical atomic MOT, 
the ground state has a larger Zeeman degeneracy than 
the excited electronic state. A quadrupole magnetic field 
gradient provides a spatially dependent energy shift of 
the G" — 1 ground state. A molecule at position 'r' 
(Fig. ID) in the upper Zeeman level preferentially scat- 
ters photons from the laser propagating to the left due 
to selection rules and laser detunings. This results in 
a restoring force towards the center of the trap. Since 
we modulate the cooling light polarization to destabi- 
lize dark states (in the G" = manifold) we must also 
modulate the direction of the magnetic field in phase 
with the light polarization to maintain a restoring force 
for the MOT. This is shown schematically in Figure ID. 
Atomic MOTs with magnetic fields oscillating at 5 kHz 
can produce stable three-dimensional trapping (22). Our 
modulation frequency of 2 MHz is set by the optical 
pumping rate. 

To fully describe the cooling and trapping forces for 
YO would involve 44 molecular levels, 15 optical fre- 
quencies with time-dependent polarization, and a time- 
dependent magnetic field. Such a calculation is beyond 
the scope of this report. Nevertheless, a simple multi- 
level rate equation model (23) can be used to extend the 
results from the two-level models to provide physical 
insight into the observed dynamics of a YO MOT. In 
the limit of small laser detunings and low laser power, 
the result from the two-level system for the Doppler and 
magneto-optical cooling force in one dimension can be 
expressed in the form (24): 

F = —f3v — nr 

-8hk 2 5s 
P ~ 7 (1 + so + {25 h) 2 ) 2 

k = fi'A/3/hk 

Here, F is the total force experienced by the molecule, 
(3 characterizes a viscous drag (Doppler) force, pro- 
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portional to the molecule's velocity, v, and k repre- 
sents a magneto-optical spring force proportional to 
the molecule's displacement, r, from the magnetic 
field minimum. The forces are parameterized by the 
wavenumber, k, of the cooling transition, the detuning, 
5, of the cooling laser from the molecular transition, the 
resonant saturation parameter s , the differential mag- 
netic moment between ground and excited state //, the 
magnetic field gradient A, and reduced Planck's con- 
stant h. 

For multi-level systems with N ground states and a 
single excited state the optimal (maximum) damping pa- 
rameter /3jv scales roughly as fi^ <~ 2fii/{N + 1), and 
is achieved when all grounds states are driven. This 
results from lower photon scattering rates due to de- 
creased population in the excited state from the addi- 
tional degeneracies of the N ground states. Despite the 
slower damping rates, the Doppler cooling limit for the 
multi-level system remains unchanged from the two- 
level result, Tdop = hj/2k-B — 116 /iK. The slower 
damping rates are balanced by slower heating rates from 
photon recoils. For our time-dependent magnetic field, 
it suffices to replace the two-level system k with a time 
averaged R. Averaging the magneto-optical force over 
a single cycle yields R = (2^/2/tt) cos((/>)k, where n is 
calculated with the root-mean-square (RMS) field gra- 
dient, and (f> is the phase between the modulations of the 
MOT field and the optical polarization. By changing the 
relative phase <j>, we can change the magnitude and sign 
of the MOT spring force independently of the Doppler 
force. 

To characterize the MOT we use a cryogenic buffer 
gas molecular beam apparatus (25), depicted in Figure 
2. The YO molecules are produced via laser ablation 
of a sintered Y2O3 pellet located inside a copper cell 
filled with a 4.5 K helium buffer gas. The YO molecules 
thermalize translationally and rotationally via collisions 
with the 4.5 K helium buffer gas. In-cell laser ab- 
sorption measurements indicate initial \X, 0; N" = 1) 
densities on the order of 10 10 cm~ 3 , corresponding to 
more than 10 10 molecules produced per ablation pulse. 
A YO molecular beam is formed by extraction of the 
molecules through a 3 mm diameter aperture in the side 
of the buffer gas cell. The molecular beam is colli- 
mated by a second aperture, 2.5 mm in diameter, placed 
130 mm from the buffer gas cell. The direction of prop- 
agation of the beam is defined as the z— axis. This re- 
sults in a molecular beam with longitudinal velocity of 
v z — 120 m/s, with full -width-half-maximum (FWHM) 



spread of 40 m/s in the longitudinal velocity distribu- 
tion (T z ~ 3.3 K), as determined by Doppler shift 
fluorescence spectroscopy. The transverse temperature 
of the molecular beam after the collimating aperture is 
Tj_ - 25 inK. 

Following the collimating aperture the molecular 
beam travels to a 10 cm long interaction region, where 
the molecules interact with the cooling lasers. The cool- 
ing lasers have a FWHM beam diameter of ~ 3 mm 
and make 1 1 round-trip passes through the interaction 
region to provide a molecule-laser interaction time of 
tint ~ 275 /US. The propagation direction of the cooling 
lasers defines the x— axis. Additionally, the magnetic 
field coil used for the MOT has a rectangular baseball 
coil geometry, with dimensions of 5 x 5 x 15 cm. The 
coil consists of 25 turns of Litz wire in series with a 
tuning capacitor, forming an LC-resonator with resonant 
frequency uj a = 2-zr x 2 MHz and quality factor Q = 77. 
Power is coupled into the MOT coils via a transformer 
designed to impedance match the MOT coil to a 50 fi 
transmission line. Assuming perfect coupling, a driv- 
ing power of 20 W corresponds to RMS field gradients 
of A = 6 gauss/cm. The MOT field coil is located in- 
side the vacuum chamber and is set in a low-outgassing, 
thermally conductive, electrically insulating epoxy and 
water cooled to avoid overheating of the coils during op- 
eration. We monitor the phase of the MOT field with a 
pick-up coil and phase lock the MOT field to the polar- 
ization modulation signal (18). 

Following the interaction region the molecules tra- 
verse a 30 cm long region for ballistic expansion. Fi- 
nally, the molecules enter a probe region where they are 
first optically pumped back into the ground vibrational 
state with a multi-pass "clean-up" beam consisting of 
v = 1,2 repump lasers. A retro-reflected probe beam, 
derived from the same laser beam used for the cooling 
transition, propagates along the 2;— axis and interrogates 
the molecules. The cycling fluorescence from the probe 
beam is collected along the y— axis and imaged onto a 
CCD camera. To extract transverse temperatures from 
the molecular beam image we fit the image profile along 
the axis to an expected functional form calculated 
from a Monte Carlo simulation of the molecular beam. 

Figure 3A shows transverse molecular beam profiles 
under various conditions. Curve (i) in Fig. 3A shows the 
unperturbed molecular beam (cooling lasers off), with a 
transverse temperature of 25 mK. If only the v" = 
cooling laser is turned on, the YO molecules are effi- 
ciently pumped into v" = 1,2 levels, leaving a depleted 
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beam signal shown as curve (ii). The curves (iii) and 
(iv) in Fig. 3A correspond to the YO molecules under 
the presence of cooling lasers detuned by S/2n = —5 
and +5 MHz, respectively. We define the detuning, 
5 = wi ascr — ujyo, as a uniform detuning of all the 
v" = cooling lasers from their respective transitions, 
with (5 = corresponding to no observable change in 
temperature. The v" = 1,2 repump lasers remain on 
resonance. Cooling of the molecular beam is observed 
as a narrowing of the molecular beam profile and in- 
crease in the number of molecules at the center of the 
beam. As expected, for small negative (positive) de- 
tunings we observe cooling (heating) of the molecular 
beam. Figure 3B shows the Doppler force dependance 
over a range of laser detunings. The observed tempera- 
tures agree well with a simulation of the force. Due to 
the closely spaced hyperfine levels in the G" = 1 mani- 
fold of the ground state, a laser with S/2tt = +10 MHz 
for the G" = 1 , F" = level will act as a negatively- 
detuned beam for the G" = 1,F" = 1,2 levels. The 
simulation of the Doppler force indicates that this effect 
is responsible for the cooling observed at large positive 
detunings. Additionally, we have negatively detuned the 
laser for each ground hyperfine manifold independently, 
while leaving the remaining v" = hyperfine lasers on 
resonance, and observed cooling in each case. This indi- 
cates that all hyperfine manifolds in the ground state par- 
ticipate in cooling. Figure 3C shows the dependance of 
the Doppler cooled temperature on polarization modu- 
lation rate for a YO molecular beam of initial transverse 
temperature of 15 mK and laser detuning S = — 5 MHz. 
The cooling and photon scatter rate is clearly limited 
by Zeeman dark states for polarization modulation rates 
less than 5 MHz. 

The final Doppler cooled beam temperatures are in 
the range of 5 to 10 mK, well above the Doppler limit of 
116 /uK. The simulations of the Doppler cooling force 
indicate this is due to a finite interaction time and cool- 
ing rate. The final temperature of the beam can be ex- 
pressed as Tf = T; x cxp[— tintTn], where 7] = 25 mK 
is the initial beam temperature, and the Doppler cool- 
ing rate T D = (2/3/m) where m = 105 AMU is 
the mass of YO. This implies an experimental value of 
Tq <~ 5 x 10 3 s _1 , which is in good agreement with 
the cooling rate To = 8xl0 3 s _1 predicted by the 
multi-level rate equation model. 

Under optimal cooling conditions we observe that 
85% of the molecules remain after cooling. This im- 
plies that the branching ratio into dark states is on the 



order of r\ <~ 10~ 4 . One possible loss channel is via de- 
cay to the A' 2 A 3 /2 electronic state, shown in Fig. 1A. 
A calculation of the dipole transition strength for the 
A 2 Tli/2 — > A' 2 A 3 / 2 transition (26) indicates that the 
branching fraction could be as large as rj = 4 x 10~ 4 . 
Molecules that decay into A' 2 A 3 / 2 will subsequently 
decay into the ground X 2 S, N" = 0,2 levels (27), 
acquiring a parity flip from emitting an additional pho- 
ton. This may be a useful method for accumulating 
molecules in the ground rotational state. Alternatively 
the X 2 S -> A' 2 A 3 / 2 transition at 690 nm could be used 
for narrow-line cooling (28) or to pump the molecules 
back into the N" = 1 level by applying a weak electric 
field to mix parity states in the A' 2 A 3 / 2 level. 

With the Doppler cooling process well characterized, 
we can now turn to the MOT's properties. Figure 4A 
shows a typical molecular beam image after passing 
through the 1-D MOT. Curve (i) shows the unperturbed 
molecular beam. Curves (ii) and (iii) show the molec- 
ular beam after passing through the 1-D MOT for rel- 
ative MOT phase = 0° (trapping) and = 180° 
(anti-trapping), respectively. While both curves (ii) and 
(iii) exhibit cooling, curve (ii) clearly shows the en- 
hancement of molecules at the center of the beam due 
to the MOT spring force. Figure 4B shows the de- 
pendance of the final temperature of the molecules on 
the relative MOT phase, <f>. The observed dependance 
agrees well with a Monte Carlo simulation of the MOT, 
yielding a value of the MOT oscillation frequency of 
Umot = <~ 2n x 155 Hz. This agrees well 

with the predicted value of wmot ^ 2tt x 160 Hz, de- 
rived from the experimentally measured value of j3 and 
calculated magnetic field gradient ^4rms = 6 gauss/cm. 
Even though the molecules traverse the magneto-optical 
trap in a fraction of a trap oscillation, the molecular 
beam intensity enhancement due to the MOT is clear 
and can be well explained by the optical pumping rate 
model. 

Figures 4C and 4D compare the observed tempera- 
tures for Doppler cooling and the MOT for varying cool- 
ing and repump powers, respectively. Simulations of the 
cooling force are shown as a solid line. In Fig. 4C we 
note that at large cooling powers, the Doppler cooling 
rate decreases (larger final temperature) due to power 
broadening of the cooling transition. Although our sim- 
ple model also predicts this effect for the MOT, we ob- 
serve the coldest MOT temperatures for the largest cool- 
ing power. This disagreement is likely due to the more 
complicated level structure in the MOT that is not ac- 
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counted for in the simulation. In Fig. 4D, at low repump 
power, the cooling rate for both Doppler cooling and the 
MOT is limited by the repump time out of the v" = 1,2 
levels, while at large repump power, the cooling rate sat- 
urates. 

While the 1 -D MOT provides simple pedagogical un- 
derstanding, to enhance the brightness of the molecu- 
lar beam we have implemented Doppler cooling and the 
MOT in 2-D. The cooling lasers then propagate along 
both the x and y directions. The action of the MOT 
along the y-direction directly resulted in an increase of 
detected molecules. The observed value of (3 X in 2-D is 
roughly half that observed in 1-D (18). This is consis- 
tent since half the photons are now scattered along the 
y-axis and do not provide any cooling in the ^-direction. 
We also note that for cooling in dimensions of 2 or more, 
modulation of the magnetic field and polarization may 
not be necessary to remix hyperfine states, simplifying 
the magnetic coil design. Type II atomic MOTs operate 
with similar level structures, and remixing occurs from 
lasers propagating in orthogonal directions (29). 

From our observations of the 1-D and 2-D MOT, 
we can estimate the properties of a 3-D MOT based 
on this system. We first define a generic 1-D damp- 
ing rate, T — (3/m. It is straightforward to show that 
the phase space density of molecules in a 3-D MOT 
will evolve as p oc nT^' 2 oc ( o exp[(3r/2 - rn P )t\, 
where r\ is the branching ratio into dark states and 7 P 
is the photon scattering rate. Phase space increase re- 
quires f] < 3r/(27 p ) ~ 10~ 3 , where we have substi- 
tuted typical parameters for YO, T = 2 x 10 3 s -1 , and 
7 P = 3 x 10 6 s _1 . The lifetime of molecules in the MOT 
will be limited by optical pumping into dark states, with 
tmot ~ (T/p) -1 > 1 ms for 77 < 4 x 10~ 4 . Typical 
capture velocities will be on the order of several m/s. 
Loading of a 3-D MOT could be achieved either from a 
slow molecular beam (14, 30) or possibly directly from 
molecules produced inside a buffer gas cell. Thus, a 
clear path exists to a 3-D MOT that produces cold, dense 
samples of diatomic molecules. 
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Figure 1: (A) YO vibronic structure. Dashed arrows indicate decay paths with corresponding Franck-Condon 
factors, q (17). Solid arrows indicate cooling and repump laser transitions. (B) Rotational and hyperfine structure 
of the X and A states. Solid arrows indicate the three hyperfine pumping components used in this work. (C) 
Zeeman structure for the X 2 Yj + , N" = 1 state. (D) Schematic of the MOT level structure and the modulation 
waveforms for optical polarization and magnetic field. 
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Figure 2: Depiction of the MOT apparatus, shown in its 1-D implementation for clarity. The 2-D system has 
cooling laser beams propagating in the vertical directions as well. YO molecules are produced via laser ablation 
(green) inside a cryogenic buffer gas cell, shown at left. The YO molecular beam (shown blue) is collimated by an 
aperture and then passes through the cooling region, shown at center. The cooling region consists of a rectangular 
magnetic field coil and a multipass laser beam interaction region. The multipass consists of a pair of mirrors and 
A/4 waveplates to provide the correct polarization of light for the MOT over the 10 cm interaction length. After 
passing through the cooling region, the YO molecules are optically pumped into the vibrational ground state (red 
beams) and the molecular beam is imaged using resonant fluorescence and a CCD camera, shown at right. Image 
credit: Brad Baxley. 
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Figure 3: 1-D Doppler cooling. (A) Molecular beam profiles for (i) an unperturbed beam, (ii) depletion due to 
lack of vibrational repumping, (iii) Doppler cooling and (iv) Doppler heating. (B) Transverse temperature vs. 
cooling laser detuning, S, for unperturbed (black) and Doppler cooled (red) beams. (C) Transverse temperature vs. 
polarization modulation rate with an initial beam temperature of 15 mK. The solid lines in (B) and (C) represent 
simulations using a multi-level rate equation model for the Doppler cooling force. 
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Figure 4: (A) Molecular beam profiles for (i) an unperturbed beam, (ii) a 1-D trapping (tfi = 0) MOT, and (iii) an 
antitrapping </> = 180° MOT. Comparisons of transverse temperature for Doppler cooled (red) and MOT cooled 
(purple) beams, varying (B) relative MOT phase, (C) relative cooling power, and (D) relative repump power. 
Simulations of the cooling forces are shown with solid lines. 
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Optical cycling transition Figure IB shows the rotational and hyperfine structure of the ground X 2 T, + state 
of YO, which is of Hund's case bps (19). The ground electronic state can be labeled by quantum numbers 
\X, v"; N" , G" , F",p), where v" indicates vibrational level and N" indicates rotational level. A strong Fermi 
contact interaction couples the electronic spin S = 1/2 to the nuclear spin I = 1/2 of the 89 Y to form an inter- 
mediate quantum number G = S + I. The electron spin-rotation interaction then couples molecular rotation N to 
form total angular momentum F = G + N. The rotational states have parity, p = (— 1) , denoted by +/— in 
Figfllb). The excited A 2 H 1 / 2 state is labeled by Hund's case a quantum numbers, and can be labeled by quantum 
numbers \A, v'; J', F',p), where J' is the total electronic angular momentum, and F = J + I. The lowest rota- 
tional state J' = 1/2 is split into two states of opposite parity separated by the A-doubling of 4.5 GHz (17). Each 
parity state consists of a pair of hyperfine manifolds F 1 = 0, 1 that are not spectroscopically resolved. 

Parity selection rules require the \ A, v'; 1/2, F', +) state to decay to a negative parity ground state, implying N" 
must be odd. Each ground level \X, v"; N" , G" , F") can be expressed as a superposition of states expressed in 
the \X, v"; J" = N" ± 1/2, F") basis. The angular momentum selection rule | J" — J'\ = 0, 1 then only allows 
decays to the \X, v" , N" = 1, G", F", -) state. Mixing of the \X, v"; 1, 1, 2) with the \X, v"; 3, 1, 2) level or the 
\A, 1/2, 1) with the \ A, 3/2, 1) level via the hyperfine interaction can lead to branching loss ratios on the order of 
(10 MHz/30 GHz) 2 < 10~ 6 and can be ignored. 

One additional complication arises when simultaneously driving the optical transitions from the energy degen- 
erate sub-levels within the hyperfine manifolds. In particular, for transitions of the the type F" — > F' = F" 
or F" — 1, where F" is an integer, optical dark states will arise for any choice of static laser polarization (21). 
In these dark states, the ground molecular hyperfine state is no longer coupled to the excited state by the opti- 
cal cooling laser, and the cooling ceases. The common solution is to destabilize the dark states by making them 
time dependent, either by applying an external magnetic field to make the energy levels non-degenerate, or by 
modulating the polarization of the driving optical field. While application of a magnetic field is straightforward 
experimentally, for this level structure in YO it does not work. The \X, v"; 1, 0, 1) state is a spin singlet state, and 
the energies of its Zeeman sublevels are insensitive to application of an external magnetic field, as shown in the 
upper panel of Fig. 1C. Therefore, we instead modulate the polarization of the cooling light between a + and a~ . 
In order to maintain optimal photon scattering rates for laser cooling, the modulation rate should be similar to the 
optical pumping rate, which in our case is on the order of several 10 6 s _1 . This is achieved experimentally with a 
voltage-controlled waveplate, also know as a Pockels cell. 

Laser cooling setup To produce the light necessary for creating a cycling transition in YO we use a setup consist- 
ing of three lasers, as shown in Figure S 1 . The laser light for the main cooling transition at A = 614 nm is generated 
from a ring cavity dye laser. The dye laser frequency is stabilized to the \X, 0; 1, 0, 1) — > \A, 0; 1/2, F' , +) tran- 
sition frequency. Cooling light for the \X, 0; 1, 1, F") manifold is generated by frequency shifting a portion of 
the main cooling beam using an acousto-optic-modulator (AOM) operated at two frequencies, 760 and 774 MHz. 
The v" = 1 and v" = 2 repump light at 648 nm and 649 nm is generated from a pair of master external cavity 
diode laser (ECDL) and injection locked slave laser setups. The 3 mW output from a master ECDL is stabilized 
to the corresponding \X, (1, 2); 1, 0, 1) — > \A, (0, 1); 1/2, F' , +) transition. This master ECDL laser is used to 
injection lock a high power laser diode (~ 80 mW). Most of the slave diode power is used for repumping of the 
\X, (1, 2); 1, 0, 1) level, while a portion (< 10 mW) is frequency shifted using a double pass AOM and used to 
injection lock a second high power diode for repumping the \X, (1, 2); 1, 1, F") manifold. To address all of the 
F" = 0,1,2 manifolds in the G" = 1 state, the AOM frequency is modulated at a rate of 1 MHz to provide a 
frequency shift between 760 MHz and 774 MHz. For each vibrational level v" = 0, 1, 2, the G" = and G" = 1 
beams are combined on separate polarizing beam splitter cubes. The v" = 1,2 repump light is then combined on 
a non-polarizing 50 — 50 beam splitter cube, with half of the power for the cooling beam and half of the power 
for the "cleanup" repump beam. The v" = 1,2 repump beam is combined with the v" = cooling beam using a 
dichroic mirror. The cooling beam at this point consists of frequencies for v" = 0,1,2 levels with the beam for 
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each level consisting of a pair of beams with orthogonal linear polarizations, for G" = and G" = 1 manifolds. 
This beam then passes through the Pockels cell allowing the polarization of the beams to be modulated between 
C7+ and o~ at rates of up to 10 MHz. 

We frequency stabilize the cooling and repump lasers via optical heterodyne measurements with light from a 
self -referenced octave-spanning erbium-doped fiber frequency comb. 

Magneto-optical trap phase stabilization During typical experimental conditions, the MOT coils are turned 
on and off, leading to cycling of the coil temperature and drifts of the resonant frequency, ojq. Since the MOT is 
driven at a fixed excitation frequency, drifts of the resonant frequency will lead to drifts in the phase shift of the 
MOT field. To stabilize this phase shift, we monitor the phase of the MOT field using a pickup coil and phase 
lock it to the polarization modulation reference signal using the circuit diagram shown in Figure S2a. Figure S2b 
shows the final MOT temperature versus the drive power of the MOT coils. The exact dependence of the final 
MOT temperature on magnetic field gradient is complicated and not fully understood. 

Comparison of 1-D and 2-D MOT Figure S3 shows a comparison of the MOT operation in one and two trans- 
verse dimensions. The total available cooling laser power is split equally between laser beams that propagate along 
the x— and y— axes. Figure S3A shows the molecular beam profile along the x— axis for 1— D cooling lasers prop- 
agating along the x— axis. Figure S3B shows the molecular beam profile along the x— axis for 2— D cooling lasers 
propagating along both the x— and y— axes. The 1-D MOT exhibits a faster cooling rate along the x— axis than the 
2-D MOT, as shown by the molecular beam's narrower width. This arises from the finite photon scattering rate. In 
two dimensions, the photons scattered by the molecules from the laser propagating along the y— axis do not lead 
to cooling along the a;— axis. They do, however, lead to cooling and compression along the y— axis, which results 
in more molecules in the field of view of the probe beam (propagating along the x— axis). Thus, the 2-D MOT 
produces a brighter molecular beam, as shown by the larger molecule number in Figure S3B. 
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Figure SI: Laser cooling setup for YO. Abbreviations: ECDL, external cavity diode laser; OI, optical isolator; 
AOM, acousto-optic modulator; A/2, half-wave retardation plate; PBS, polarizing beam splitter. Due to the limited 
output power of the ECDL, we use injection seeding to boost the usable optical power for the MOT experiment. 
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Figure S2: (A) Circuit diagram for the MOT coil. The phase of the MOT is monitored via a pickup coil. The signal 
from the pickup coil is mixed with the reference signal from the polarization modulation and used to phase lock 
the MOT coil to the polarization modulation. (B) The measured molecular MOT temperature as a function of RF 
drive power, which corresponds to the MOT magnetic field gradient under modulation. 
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Figure S3: Comparison of 1-D and 2-D MOTs. (A) Molecular beam profiles for an (i) unperturbed beam and (ii) 
1-D MOT. (B) Molecular beam profiles for an (i) unperturbed beam and (ii) 2-D MOT. The 1-D MOT exhibits a 
faster cooling rate along the observed axis, while the 2-D MOT produces a brighter molecular beam. 
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